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We carry out a feasibility study toward a lattice QCD calculation of the neutron electric dipole mo- 
ment (NEDM) in the presence of the 6 term using two different approaches. In the first method, 
we calculate the CP-odd electromagnetic form factor F3, which becomes the NEDM in the zero 



momentum transfer limit. At the first order in 0, we derive a formula connecting the lattice 
three-point function to the CP-odd electromagnetic form factor. In the second method we directly 
extract the NEDM from the energy difference between spin-up and spin-down neutron states in 
the presence of a constant electric field, without expanding a small but non-zero 6. We test both 
approaches numerically, employing the domain-wall quark action with the RG improved gauge 
action in quenched QCD at a -1 ~ 2 GeV on a 16 3 x 32 x 16 lattice, and further applying the 
second method to the clover quark action at a similar lattice spacing and nucleon mass. We ob- 
tain good signals from both approaches. In particular the second method works well with both 
fermion formulations. 

XXIIIrd International Symposium on Lattice Field Theory 
25-30 July 2005 
Trinity College, Dublin, Ireland 



* Speaker. 

' shintani@het.ph.tsukuba.ac.jp 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos . sissa . it / 



Neutron electric dipole moment on the lattice 



Eigo Shintani 



1. Introduction 

One of the most stringent constraints on possible violation of parity and time-reversal sym- 
metry in the strong interaction comes from the measurement of the electric dipole moment (EDM) 
for neutron (NEDM) d n and proton (PEDM) d p . The current upper bound is given by \d n \ < 
6.3 x 10" 26 e-cm(90% C.L.)[§] and \d p \ < 5.4 x lfr 24 e-cm[|J. On the other hand, QCD allows 
a gauge invariant renormalizable CP-odd term in the Lagrangian. Crude model estimations [|], ||] 
and the experimental bound on \d n \ yield a very stringent bound B < 0(1O~ 10 ). 

Current model estimates give not only an order of magnitude different results on d n but also 
differ even in its sign. Clearly a first principles determination of d n from lattice QCD is required 
to determine the value of B, if a non-zero value is found for NEDM in future experiments. Indeed 
serious attempts toward lattice QCD calculations of NEDM have just been started recently [JI], 0, . 
In this report we present our results on NEDM calculations based on two different methods, one of 
which has already been published^]. 

2. NEDM from the electromagnetic form factor 
2.1 Formulation 

Let us consider the electromagnetic form factor of the nucleon defined by 

~Hq 2 



(N(p,s)\J™\N{p',s'))=u(p,s) 



-q v o^ v Y 5 -\ 



u(p>,s'), (2.1) 



2mn 

where /™ is the electromagnetic current, q = p — p' is the momentum transfer, \N(p,s)} is the 



on-shell nucleon state with momentum p, energy po = y m$j + p 2 and helicity s. On the right- 
hand side we explicitly write down the form factor Fj only, which is related to NEDM as \d n \ = 

lim F$(q )/2mN = F^(0)/2mN . 

q 2 ^0 

The electromagnetic form factor can be extracted from the following three-point correlation 
function: 

G e NhN {q,t,%) = e (N(p,t)J* M (q,z)N(p',0)) e 

= e-V( ^ - T )e- £ ^ ^ £4(p, s )4(p,^(^)4(p^ S ')4(/J , ,/) + ■ • • , (2.2) 

<(<?) = y^{q 2 ) + ^^-q v O^ + id t^lq v y 5 a^ v + F A (q 2 )(q^- q^Ys) + 0(B 2 ) 

(2.3) 

where N(p,t) and N(p,t) are the interpolating fields of the nucleon at time t, and F\ and F2 are 
CP-even form factors while F 3 and F A are CP-odd. The on-shell spinor satisfies 



l^u N (p,s)u N (p,s)=Z N — =Z N — + iBZ N ——y 5 (2.4) 



and for small B we have an expansion of the forms m N e = + 0(B 2 ), 2^ = 2^ + 0(B 2 ) and 
f N (B)=f^B + 0(B 3 ). 



mil 



Neutron electric dipole moment on the lattice 



Eigo Shintani 



In lattice simulations, the three-point function is calculated as 



G e NJtiN (q,t,z) = J ^U^Y^We SQCD+idQ = G NJtlN (q,t,z)+idG^ N (q,t,T) + 0(d 2 ).(2.5) 



Comparing eq.( [2.2[ ) with eq.( |2.5[ ), we obtain two independent formulae which allow an extraction 
of F3. They are given by 



q 



2E N m N 



Hq 2 ) + 



\Z N \ 2 e- EN ^e- E '»' 

E N + m N 2 

F\{q ) + 



tr pNj iN {qJ^)iTAy5yi\ 



2E N 4m N E, 



\Z N \ 2 e- E ^-^e- E '^ 



-F 2 (q 2 



E N + m N 
2EnMn 



qi -F { (q 2 )- qi{E A N + 3mN) F 2 (q 2 ] 



fh 



(2.6) 



(2.7) 



2En Am^E^ 

where r4 = (1 + 74) /2. The coefficient fy, which we need to subtract off the mixing terms in the 
above formulae, can be extracted from the two-point correlation function as 

r >2„-E N tfN m N r 



(N(p,s)N(p,s)Q) 



2E, 



N 



"75, 



(2.8) 



and Fi and F2 can be obtained from Gnj^n- 



2.2 Numerical result 



We calculate two- and three-point functions on 730 quenched gauge configurations generated 
with the RG-improved gauge action at j8 = 2.6, corresponding to aT x = 1.81(4) GeV determined 
from p meson mass, on a 16 3 x 32 lattice. For the quark action, we employ the domain-wall fermion 
with the fifth dimension length N s = 16 and the domain-wall height M = 1.8, taking the quark mass 
of mj = 0.03 corresponding to nips/my = 0.63. The topological charge Q is calculated with the 
cooling method with an 0{a 2 ) improved definition. As shown in ref.[|I|], 730 configurations are 
sufficient to extract the CP-odd phase parameter from the two-point correlation function, and 
we obtain f N = -0.247(17). 

Subtracting the contribution composed of F\ , Fi_ , from the three-point correlation function, 
we obtain two independent estimates for F3 according to eqs^Zjj) and (2J_), which are plotted 
in Figjl] as a function of t for x = 6 being fixed. The two estimates agree at large enough t, 
demonstrating the correctness of our formulae. After fitting the results from eq^I^) at 12 < t < 15, 
we obtain 

' -0.024(5) e-fm neutron 



F 3 {q l ~0.58GeV z )/2m,v 



+0.021(6) e-fm proton 



(2.9) 



3. NEDM from the energy shift due to the electric field 

In the previous section, using the quenched lattice QCD calculation, we successfully obtained 
the signal for the CP-odd form factor F3. However, calculating NEDM still requires the limit 
q 2 — > 0, which is not so easy in lattice QCD since momenta can be varied only discretely in a finite 
box and the statistical error increases as the momentum increases. Therefore, we have investigated 
an alternative method proposed in ref.|Q], which does not require the q 2 — > extrapolation. 
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Figure 1: The F3 term for neutron (left) and proton (right) as a function of t with the current inserted 
at T = 6. Circles represent the result with 1^475 projection in eq.( |2.6[ ), while triangles that with T^jj in 



eq.(2.7), averaged over j = 1,2,3. Straight lines show the fitting results for the former. 



3.1 Formulation 

A non-zero NEDM generates an energy shift of the nucleon in the presence of a uniform and 
static electric field E. For E = (0,0, E), for example, the energy difference between the spin-up 
state iVf and the spin-down state iVj takes the form m N e (E) — m N e (E) = d®E + 0(E 3 ) . 

We include a uniform electric field by the substitution Us(x) — » U?,(x)e qEt , where the quark 
charge q is 2/3 for the up quark or — 1 /3 for the down quark. This substitution, necessary for a real 
electric field in Minkowski space, violates the periodicity in the t-direction. In order to minimize its 
effect at the boundary, we have to take the electric field as small as possible. Choosing appropriate 
spin components in the neutron propagator, we have 



R(E,6,t) = l(N 1 N 1 )e(E,t)}/[(N 2 N 2 )e(E,t)]=Ze- d ° Et+0 ^ + ---, 



where 



(N a N a )e(E,t) = (N a (t)N a (0)e i6Q ) E 



(3.1) 



(3.2) 



and \&)e represents the vacuum expectation value of 6 without the 6 term but with E. Since the 
neutron propagator calculated at 6 = is averaged with the weight e' e ® ("reweighting") to obtain 
the 0-average, a good sampling of the topological charge Q is crucially important for a success of 
this method. 

3.2 Numerical results with the domain-wall fermion 

We take the same simulation parameters as in the previous section except for a heavier quark 
mass, nif = 0.12, corresponding to c± 2.2 GeV. The number of configurations is increased to 
1000. 

On the left of Fig.| the ratio eq.(fO]) is plotted as a function of t for E = ±0.004 and 8 = 0.1 
(circles), and for E = 0.004 and 6=0 (bursts). The behavior of the ratio is consistent with the 
expectation that the exponent depends on the sign of E and that the ratio is unity for 6 = 0. On the 
right in Fig.|| we show the effective mass of the ratio, from which we determine the fitting range to 
be 5 < t < 9. 
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Figure 2: (Left) The ratio eq.(|3.1|) as a function of t foiE = ±0.004 and B = 0.1 (circles), and for£ = 0.004 
and 0=0 (bursts). (Right) An effective mass plot of the ratio. 

As is seen on the left of Fig.||, the exponent of R(E, d,t) has a linear behavior in E as expected. 
Finally, in right of Fig.[| df r is plotted as a function of 6. A slope obtained by a linear fit becomes 



d d n /G = -0.0183(60) e • fin, d? /d = 0.0176(72) e • fm. 



(3.3) 



These values are smaller than those obtained with the CP-odd form factor at q 2 ~ 0.58 GeV 2 ( 



eq.([2.9|) ). Note, however, that the quark mass is a factor four larger in the present calculation. 



3.3 Preliminary results with the clover fermion 

Since the Wilson-clover fermion action requires much less computational costs, and there 
exists a large number of full QCD configurations generated with this action, it is important to 
test if the method of energy difference works with it. We calculate the ratio R(E,6,t) using the 
gauge configurations generated with the same parameters, and choosing the hopping parameter 
K = 0.1320 which yields a nucleon mass similar to the domain- wall case. In Fig.|] we present 
our preliminary result for the effective mass. It looks very similar to the one with the domain-wall 
action. Our results suggest that the difference of the two fermion actions does not affect the success 
of this method. 
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Figure 3: (Left) E dependence the ratio R(t,E, 6). (Right) 6 dependence of NEDM over lattice spacing a. 
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Figure 4: An effective mass plot as Fig. 2 with the clover fermion at K = 0.1320, 9 = 0.1 and E = 0.004. 

4. Conclusion and outlook 

In this report we have proposed and tested two different approaches for the calculation of 
NEDM on the lattice. In the first method the CP-odd electromagnetic form factor is employed, 
while in the second one the energy shift is calculated in the presence of an electric field. We 
have shown that both methods work well with the domain-wall quark action. Furthermore we 
have demonstrated that the second method also works with the clover quark action. We are now 
calculating NEDM at lighter quark masses with the clover quark action in quenched QCD, in order 
to see how errors increase as the quark mass decreases. Our final target is the calculation of NEDM 
on full QCD configurations generated with the clover quark action by the CP-PACS collaboration. 

This work is supported in part by Grant-in-Aid of the Ministry of Education (Nos. 13640260, 
15204015, 15540251, 15740134, 15740165, 16028201, 16540228, 17540249, 17540259 ). 
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